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Abstract. Extraction of the nucleon's strange form factors from experimental data requires a quantita- 
tive understanding of the unavoidable contamination from isospin violation. A number of authors have 
addressed this issue during the past decade, and their work is reviewed here. The predictions from early 
models are largely consistent with recent results that rely as much as possible on input from QCD sym- 
metries and related experimental data. The resulting bounds on isospin violation are sufficiently precise 
to be of value to on-going experimental and theoretical studies of the nucleon's strange form factors. 

PACS. 13.40.Gp Electromagnetic form factors - 14.20.Dh Protons and neutrons - 12.39.Fe Chiral La- 
grangians - 12.39.Ki Relativistic quark model 



1 Motivation 

Isospin violation is generally a small effect. For example, 
consider the nucleon mass splitting, (m n — m p )/m p = 
0.1%. One similarly expects isospin violation to have a 
small impact on the nucleon's electromagnetic and weak 
form factors. However, this does not imply that isospin 
violation must be small relative to strangeness effects. To 
illustrate, recall that an explicit calculation in the elec- 
troweak theory leads to 

G P x Z (q 2 ) = 

(1 - 4sin 2 e w )G p x (q 2 ) - G x {q 2 ) G x (q 2 ) G u /(q 2 ) (1) 

for electric (X = E) and magnetic (X = M) form factors. 
Experimental studies 1,2,34,5,6 show that the sum of 
the last two terms on the right-hand side is small. The 
size of isospin violation, G x ^'{q 2 ), relative to strangeness, 
G s x (q 2 ), is not obtained from these experiments. 

In what follows, theoretical studies of G u x d {q 2 ) will be 
reviewed [7|l8]l9|ll()llll| . (Our entire discussion of isospin vi- 
olation also fits within the more restrictive category called 
"charge symmetry breaking" and that language is used, 
for example, in Ref. [B].) If the current understanding of 
these isospin violating effects is sufficiently precise, then 
the data from Refs. 1,2,3,4,5,6 allow for a determination 
of the authentic strange quark effects, G s x (q 2 ), which are 
of great interest to many people at present. 

Independent of any chosen theoretical approach, each 
isospin violating form factor is simply the difference of 
isoscalar and isovector terms, 

G u x d {q 2 )^G x (q 2 )-G x (q 2 ), (2) 

where G x is obtained from 

(p |u7 M w + d7 M e?| p) - (n |u7 M u + J-7 M rf| n) (3) 



and G x is obtained from 

(p |m7 m m - d7 M d| p) + (n {uj^u - dj^dl n) (4) 

in a straightforward manner (see Refs. [7llll| for details). 
Furthermore, we know that all isospin violation is ulti- 
mately a consequence of unequal quark masses, m u ^ mj, 
("strong breaking") and unequal quark electric charges, 
e u 7^ ed ("electromagnetic breaking"). The task for each 
theoretical approach is to determine the combinations of 
nucleon matrix elements shown in Eqs. J2J and Q, with 
both types of breaking included. 

Since the sum of strangeness and isospin violation in 
Eq. is measured to be a small fraction of the total form 
factors, and since isospin violation itself is expected to be 
a small fraction of the total form factors, it is reasonable 
to neglect contributions containing both strangeness and 
isospin violation as doubly (i.e. negligibly) small. This al- 
lows G^ d (q 2 ) to be calculated without dynamical strange- 
quark effects. Such an approach is clearly advantageous for 
chiral perturbation theory, where addition of a dynamical 
strange quark leads to severe degradation of convergence 
properties of the chiral expansion. All of the theoretical 
studies to date have computed isospin violation without 
dynamical strange quarks. 

Our discussion will be approximately chronological. 
The constituent quark model studies of Dmitrasinovic and 
Pollockjjj and MillerjH], are discussed in Sec. [21 followed 
by Ma's use of a light-cone meson-baryon fluctuation mod- 
el [HI in Sec. [3 The constraints of chiral symmetry are dis- 
cussed in Sec. 01 based on a collaboration with MobedfTTi] 
using chiral perturbation theory. Section [SI reviews the re- 
cent results from work with KubisjTT] that combines chiral 
perturbation theory with resonance saturation and infor- 
mation from dispersion analyses. The final section, Sec. [HI 
provides a brief summary. 
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2 Constituent quark model 

In a constituent quark model, we might expect the scale 
of strong breaking to be set by constituent quark masses, 
which the authors of Ref . take to be (mjj — mu ) / rriQ ~ 
1.2%, and the scale of electromagnetic breaking to be set 
by a ~ 0.7%. To do better than this order-of-magnitude 
guess, an explicit calculation is required. 

The first calculation was carried out by Dmitrasinovic 
and Pollock^ using oscillator confinement and a Coulomb 
potential, 



H — Hq + Vem 
3 



^ 2m< 2 ^ v °' ' 



i=l 
3 



i<j 



4tt |ri — rj| 



(5) 
(6) 

(7) 



Recall that this choice for Ho produces Gaussian spatial 
wave functions. With parameters fixed to be mq — 330 
MeV, mo— mu = 4 MeV, and k determined from the ex- 
perimental mass difference between Z\(1232) and nucleon, 
Ref. |Zj finds 



<4(0) = <4(0) 

(4) 

(r 2 M ) 



m D 
mu 



m v 
m D 



<(0)-G^(0) 
1.1%, 

0.4% , 



0.008 , (8) 
= 0, (9) 
(10) 

(11) 



where r\ denotes a squared radius as usual. These esti- 
mates are compatible with our order-of-magnitude guesses, 
but what confidence level should be assigned to the precise 
values? 

Here is a list of some limitations of this model (all of 
which are mentioned explicitly in Ref. [7]): 

— The chosen parameters lead to, m n — m p = 3 MeV, 
230% above experiment. 

— The nucleon charge radius, Vvb) = 0-62 frn, is 30% 
below experiment. 

— Gaussian spatial wave functions cause both Ge{<1 2 ) 
and Gm{<1 2 ) to be unrealistic at large q 2 . 

— The strong hyperfine interaction is omitted. 

— Chiral symmetry is absent. 

In Section IV B of Ref. 0, the authors conclude, ". . . we 
may expect to have calculated the correct sign and order 
of magnitude of the effects of interest." 

Subsequently, Miller[H] used a more complete constitu- 
ent quark model and explicitly addressed each of the lim- 
itations listed above, except chiral symmetry. (Miller says 
chiral symmetry is implicit in the charge symmetry con- 
serving pion cloud of this model^Jl-) His Hamiltonian is 



H = K + V cl 



9 i 



(12) 



Table 1. Parameter values for the three Gaussian models of 
Ref. 0. 



Model 


1 


2 


3 


VP (fin) 


0.7 


0.6 


0.5 


OLs 


2.3 


1.2 


0.35 


Tit o — mu (MeV) 


5.2 


3.8 


2.3 



1=1 



Pi 

2m, 



(13) 



i;, H <>^J </.'/, ( — ^ ) 

i<3 



2 | 4 ff(i). ff Q-) 

m 2 3 m 2 



(14) 



i<3 



1 1 , 4a(i)-a(j) 



mf m~ 3 TYliTYlj 



(15) 

Two options for V^on are studied; one gives Gaussian form 
factors, 



W(p, A) = Nexp 



1 



P 2 + A ; 
-2/3 



P=^(rr-r 2 ), 

A= -^=(r 1 +r 2 -2r 3 ), 
and the other gives power-law form factors, 



(16) 
(17) 
(18) 



G E (Q 2 ) = 



A 2 



) 2 + A 2 



(19) 



For the Gaussian case, Ref. [Hj uses the experimental 
proton magnetic moment to determine m — 337 MeV, 
the experimental tra — m^ value to fix a s as a function 
of (3, and then considers three choices for (3 which define 
three models as shown in Tabled Figures 4, 5, 6 and 7 of 
Ref. |] display the resulting effects of isospin violation on 
Ge(q 2 ) and Gm(q 2 ), including the separate contributions 
from K, V cln and V g . In each case, the isospin violating 
contributions vanish for q 2 = 0, and are less than 0.2% 
in magnitude for a momentum transfer of 0.1 GeV 2 . For 
power-law form factors, Ref. |Hj obtains A — 5.90/fm from 
the experimental value of (r E ), and all isospin violating 
effects remain small. 

The combined work of Refs. [7||S] provides an excel- 
lent understanding of isospin violation within the con- 
stituent quark model, but we must now ask which fea- 
tures of the results are a true reflection of nature, and 
which are model-dependent. For example, there are sym- 
metries in these quark models that lead to G^f(0) = 0, 
recall Eq. but, as we'll see in Sec. 01 this is not a sym- 
metry of nature. It is therefore interesting to explore other 
theoretical approaches as well. 
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3 Light-cone meson-baryon fluctuation model to arrive at 



To discuss the method used by Ma|§], we must think 
at the hadron level rather than the constituent quark 
level. Fluctuations of a nucleon into a virtual pion-nucleon 
pair, p — > n + n and n — > ir~p, occur commonly. As al- 
ways, isospin violation arises from strong breaking (at the 
hadron level, this means m„ — m p ^ 0) and from elec- 
tromagnetic breaking (which now means the Coulomb at- 
traction of a charged/charged ir~p pair is different from 
the charged/neutral 7r + n pair). The expression for isospin 
violation in the magnetic form factor at vanishing momen- 
tum transfer, is 



GT(0) 



P(n - ix-p) - P(p - Tr+n)) (V p 



(20) 

where P() denotes a fluctuation probability, and is 
the magnetic moment for the neutron's fluctuation. 

To determine begin with the fact that total an- 

gular momentum is orbital angular momentum plus pro- 
ton spin, 



1 1 

2' 2 



1,1 



This leads directly to 



3 3 
4.75 



---) - 

2' 2/ s 



2m n 
m n - + nip 



1,0 



1 1 

2' 2 



3 

21) 



fJ-n 



-2.84. (22) 



To determine the difference of fluctuation probabilities 
needed for Eq. H20|) . a light-cone Gaussian wave function 
is used, 

xP = Ae -{M--ml)/(^)^ (23) 

where 



M 2 = E 



(24) 



is the invariant mass of the meson-baryon state. (Note 
the implicit assumption that A is independent of which 
nucleon is fluctuating; Ref. [5] points out that this as- 
sumption is not required, but the uncertainties associated 
with relaxing it are difficult to estimate.) From the exper- 
imental Gottfried sum rule, one finds 



P(p -> 7T+n) « P(n -> ir~p) w 0.15 . 



(25) 



The final remaining parameter is the radius, a, and it 
leads to a large uncertainty. Ref. [Hj uses two bounds, 



a = 300MeV 



and 



P{n = n~n) - P(p = 7r + ri) = 0.2% , 

(26) 



a(n -> n~p) = 205MeVl 
a(p -> n+n) = 200MeV J 

=>■ P{n = Ti-n)- P{p = TT+n) = 3%, 

(27) 



G^(0) « 0.006 -> 0.088. 



(28) 

However, Miller suggests that (205-200)/200=2.5% is too 
large for a Coulomb effect, and states: "A reasonable es- 
timate of the effect could be 0.03 nuclear magnetons." [Hj- 
Thus, the light-cone meson-baryon fluctuation model 
provides an explicit example of a model that does not 
lead to G^f(Q) = 0. In particular, this model predicts 
G U M d (0) > 0. 



4 Chiral perturbation theory 

Chiral perturbation theory is not a model; it merely ad- 
ministrates the global symmetries of QCD. This has the 
advantage of retaining no model dependence in predic- 
tions, and it has the disadvantage that quantities not 
determined by global symmetries remain completely un- 
known. In Ref. |10j . isospin violation in the nucleon's vec- 
tor form factors is studied within chiral perturbation the- 
ory to leading order for the electric case, and to next-to- 
leading order for the magnetic case. 

Because isospin violation arises from both strong break- 
ing and electromagnetic breaking, our chiral perturbation 
theory must include dynamical photons as well as dynam- 
ical pions. A single baryon will flow through the process, 
though its identity may change from proton to neutron 
or even to Z\(1232) by absorbing or emitting a pion (re- 
call the fluctuations discussed in Sec.|3|). Ref. ^Oj uses the 
formalism of heavy baryon chiral perturbation theory [Hfl 
114) . and the perturbative expansion is defined in powers 
of electric charge e, momentum q/A, pion mass m^/A, 
and mass difference (mj — mjv)/ A, where A denotes ei- 
ther AttFtt or tun- Writing the Lagrangian in standard 
notation, 



C 



ChPT 



£(D + £(2) + £ (3) + £ (4) + £ (5) + 



(29) 



Ref. 10 shows that the isospin violating form factors be- 
gin in £( 4 ) (plus the corresponding loop diagrams), and 
next-to-leading effects are in (plus the correspond- 
ing loop diagrams). For a general observable at this high 
order, two-loop diagrams routinely appear and the most 
general effective Lagrangian contains literally hundreds 
of low-energy constants (i.e. parameters whose numeri- 
cal values are not constrained by global symmetries, and 
hence unknown to chiral perturbation theory). 

For the specific case of G u ^{q 2 ), Ref. ^H] shows that 
the situation is much simpler than for a general observable 
due to five key observations: 

— G^f(q 2 ) receives no two-loop contributions up to next- 
to-leading order. (In particular, "photon+photon" loop 
diagrams lack the required spin operators, "photon+ 
pion" loop diagrams sum to anticommutators of Pauli- 
Lubanski operators, and "pion+pion" loop diagrams 
have only higher-order isospin violation.) 

— All one-loop photon effects in (q 2 ) can be absorbed 
into the physical value of m n — m p . 
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— Explicit Z\(1232) effects are found, numerically, to be 
smaller than systematic uncertainties. 

— G^f(q 2 ) contains just one combination of low-energy 
constants at leading order — it's a simple additive 
constant — and no additional low-energy constants at 
next-to-leading order. (The fact that no global sym- 
metry forces this constant to vanish is not realized in 
some constituent quark models, as discussed in Sec. [3) 

— All isospin breaking in loops, up to next-to-leading or- 
der, is ultimately due to m n — m p . 

Based on these observations, a parameter-free prediction 
for the difference G u ^{q 2 )— G^f(0) is obtained in Ref. |10| . 

Isospin violation in the electric form factor is discussed 
in Ref. ^5], but only to leading order because in this case 
there are two-loop contributions and additional unknown 
constants at next-to-leading order. 

Though it is useful to have these "pion cloud con- 
tributions" computed within chiral perturbation theory, 
the missing combination of low-energy constants indicates 
that chiral perturbation theory alone will not meet all of 
our goals, and we must turn our attention to a method 
for understanding the physics that resides within the low- 
energy constants. 



■N 



0) 



■N 



low-energy constant 



resonance exchange 



Fig. 1. The basic idea of resonance saturation for a nucleon's 
vector form factor. 



CO 



CO 



Fig. 2. Resonance saturation contributions to isospin violation 
in the nucleon's vector form factors. 



5 Chiral perturbation theory with resonance 
saturation 

The chiral perturbation theory calculations of Ref. [THj 
were reproduced in Ref. using two separate formalisms: 
a repeat of the heavy baryon calculation, and a newer 
method known as infrared regularization|T?)] . This newer 
approach is simpler to manage and needs fewer Feynman 
diagrams, but physical results must be identical. Ref. 
finds and corrects some errant factors of 2 in Ref. , and 
confirms all of the features mentioned in Sec. 21 

With the results of chiral perturbation theory now 
firmly in hand, Ref. addresses the issue of unknown 
low-energy constants. They are not constrained by chiral 
symmetry, nor by any other global symmetry of QCD. 
Since chiral perturbation theory is the low-energy effec- 
tive theory of QCD, where pions (and in our case, pho- 
tons too) are dynamical, the low-energy constants must 
account for all of the higher-energy QCD dynamics. Most 
important would be the exchange of the lightest reso- 
nances that are too heavy to appear explicitly in the ef- 
fective theory: for vector form factors, these are the p and 
ijj. Figure sketches this notion of resonance saturation. 
Mathematically, it is easy to see that the exchange of a 
heavy particle is a low-energy constant at leading order, 



propagator < 



M 2 



1 

If 2 



O 



q 

M 2 



(30) 



Of course these resonances have to be coupled in a way 
consistent with chiral symmetry. 

This type of resonance saturation was shown long ago 
to work very well in the meson chiral Lagrangian |17U181 
HHj . though it has been shown that the nucleon's vector 



form factors require inclusion of extra resonances beyond 
merely the p and u> 1201 • Nevertheless, our present goal is 
only to compute the isospin violating pieces of these form 
factors; they come from p — oj mixing, and Ref. 1 1 1| points 
out that any effects of higher resonances are more severely 
suppressed in this case. 

The contributions of resonance saturation to the isospin 
violating form factors are shown diagrammatically in Fig. [21 
The corresponding equations are 



SG u /(q 2 



5G v ' d '" 2 



0„ 



M V (M 2 



1 



ci 2 

9uF p q 



- 1 



k p M 2 
Am 2 N 



M v {M 2 -q 2 ) 2 
- (q 2 + k p M 2 ) g p F L 



9 P Fuq 2 

[(q 2 + K^Ml) 9uJ F p 



(31) 



(32) 



where My, F p and F u are the vector meson mass and de- 
cay constants, and the p — lo mixing parameter has been 
determined from experimental masses and branching ra- 
tios to be|21j 



O pu] = (-3.75 ±0.36) x 10 _3 GeV - 



(33) 



If we had numerical values for the couplings of vec- 
tor mesons to nucleons, g p , k p , and k^, then Eq. (|32|1 
could be expanded in powers of q 2 to obtain the desired 
low-energy constant (at q 2 = 0) as well as additional con- 
tributions that are technically of higher order in the chiral 
expansion. We immediately see that Eq. I|31[) vanishes at 
q 2 — as required by global symmetries, but nonzero 
contributions from higher orders in the chiral expansion 
would be obtained. 
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0.04- 



0.03- 




Fig. 3. Isospin violation in the electromagnetic form factors, 
as obtained from Ref. 



with minimal model-dependence by using dispersive anal- 
yses. The results in Fig. [3] represent a conservative deter- 
mination of isospin violation, as obtained from worst-case 
error bars. 

Because it is the sum G s x (q 2 )+G u ^_(q 2 ) from Eq. Q 
that is measured in experiments, Fig. |21 provides a theo- 
retical error bar for the extraction of G s x (q 2 ). As shown 
explicitly in Table III of Ref. jllj , modern experiments are 
already approaching this level of precision. 
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In principle, there is a wide assortment of techniques 
for attempting to quantify g pi n p , g u and k u , but finding 
an approach with rigorously-quantifiable uncertainties is 
difficult. In an attempt to rely as directly as possible on 
experimental data rather than models, Ref. jn] uses val- 
ues extracted from dispersive analyses of nucleon electro- 
magnetic form factors. For the p couplings, these analyses 
must account for the non-resonant two-pion continuum 
in addition to the Breit-Wigner p resonance. Data from 
Refs. 22,23 lead to the ranges 

4.0 < g p < 6.2 , (34) 
5.1<k p <6.8. (35) 

For the u couplings, it is sufficient to consider only pure 
zero-width resonance pole residues, and Refs. |2~2l24ll25| 
produce the ranges 

41.8 < g u < 43.0, (36) 
-0.16 < k u < 0.57. (37) 

The final results of Ref. JT], reproduced in Fig.|3 include 
the worst-case errors bars obtained from spanning the 
ranges in Eqs. (|34I37|I above; the poorly-known n u dom- 
inates the uncertainties. Notice that the range of G^f(Q) 

does not include zero, while G^' d (0) = is required. Both 
form factors are positive over the momentum range con- 
sidered. 

6 Summary 

Early studies of isospin violation in the nucleon's vector 
form factors led to a clear understanding within specific 
quark models. Use of chiral perturbation theory avoids 
all model dependence, but leaves some parameters unde- 
termined. Phenomenologically, those parameters are sat- 
urated by resonances, and numerical values are obtained 
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